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Abstract 



We show that the effective decay rate of Zeeman coherence, generated in 
a 87 Rb vapor by linearly polarized laser light, increases significantly with the 
atomic density. We explain this phenomenon as the result of radiation trap- 
ping. Our study shows that radiation trapping must be taken into account to 
fully understand many electromagnetically induced transparency experiments 
with optically thick media. 
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Coherent population trapping (CPT) was first observed in experiments establishing Zee- 
man coherence in sodium atoms . In these experiments, explained in terms of a three-level 
A-type level scheme, a laser field was used to create superpositions of the ground state sub- 
levels. One of these superpositions, referred to as the "bright" state, can interact with the 
laser field while the other superposition does not and is referred to as the "dark" state 0. 
All the population in the system is eventually optically pumped into the dark state, and 
resonant absorption of the electromagnetic field almost disappears. This phenomenon is one 
manifestation of electromagnetically induced transparency (EIT) ||.||]. 

EIT is particularly interesting because it offers a wide variety of applications ranging 
from lasers without population inversion |4|,U to new trends in nonlinear optics 0. One 
of the most striking phenomena connected with EIT is that atoms prepared in a coherent 
superposition of states can produce a steep dispersion and a large resonant index of refraction 
with vanishing absorption . Preparation of matter in such a state (which has been dubbed 
"phaseonium" ||) provides us with a new type of optical material of interest both in its own 
right, and in many applications to fundamental and applied physics. 

A common condition for applications of EIT is a high optical density of the resonant 
medium. For example, in experiments demonstrating enhancement of index of refraction 
(x' ~ 10~ 4 ) the density of the particles was N ~ 10 12 cm -3 Q. Also, for highly sensitive 
magnetometry based on atomic phase coherence the density of atoms is estimated to be 
JV~5 x 10 12 cm- 3 f§. 

For optically thick media, reabsorption of spontaneously emitted photons can become 
important. This process, called radiation trapping, has been studied extensively in astro- 
physics, plasma physics, and atomic spectroscopy || . Radiation trapping has been predicted 
and demonstrated to have a destructive effect on the orientation produced by optical pump- 



ing pOH12| . Because the spontaneously emitted photons are dephased and depolarized with 
respect to the coherent fields creating the atomic polarization, the effect of radiation trapping 
can be described as an external incoherent pumping of the atomic transitions JT2J . Under 



the conditions of EIT, there are not many atoms undergoing spontaneous emission. How- 



ever, these spontaneous photons destroy the atomic coherence in the same way as incoherent 
pumping. This effect can change the results of CPT and EIT experiments significantly 

In this Letter we report the observation and analysis of the increase of the effective decay 
rate of Zeeman coherence due to radiation trapping. We create the coherence between ground 
state sublevels of the Dx line of 87 Rb vapor, and study the effect of CPT as a function of the 
vapor density using polarization spectroscopy based on nonlinear magneto-optic rotation 
(NMOR) fp~3| P~5p . We see that the the relaxation time of the coherent state is determined 
not only by the time-of-flight of the atom through the laser beam, but also by the density 
of the atomic vapor. For atomic densities Af^5x 10 12 cm -3 the effective coherence decay 
rate increases by several times compared with the decay rate for JV » 5 x 10 11 cm -3 . This 
demonstrates the importance of radiation trapping for experiments with optically thick 
coherent media, and the need to account for it to understand the experimental results. 

Our experimental setup is shown schematically in Fig. |l| An external cavity diode laser 
is tuned in the vicinity of the F = 2 — > F' = 1 transition of the 87 Rb Dx line. The laser 
beam passes through a high-quality polarizer P 1; and then through a cylindrical glass cell 
containing isotopically enhanced 87 Rb. The length of the cell is L = 5.0 cm and its and 
diameter is D = 2.5 cm. The laser power after the polarizer is P = 2.5 mW, and the beam 
diameter is d = 2 mm. The glass cell is placed inside a two-layer magnetic shield to suppress 
the laboratory magnetic field. A homogeneous longitudinal magnetic field is created by a 
solenoid mounted inside the inner magnetic shield. The density of the Rb vapor is controlled 
with the temperature of the cell. A second polarizer P 2 (a polarizing beam splitter) is placed 
after the cell and is tilted at 45° with respect to the first polarizer. The two beams emerging 
from the polarizing beam splitter are detected with detectors D\ and D2. A simple analysis 
of the signals from the two channels gives the angle of rotation of the polarization <p and the 
transmitted intensity J out . 

We consider the linearly polarized light as two circular components E + and £L which 
generate a coherent superposition of the Zeeman sublevels (a dark state). To study this 
dark state we apply a longitudinal magnetic field in the direction of light propagation which 
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leads to a splitting of the \b±) states of magnitude Mo/2 = —jisB where B is the mag- 
netic field, hb is the Bohr magneton. Because the two circular components interact with 
ground-state sublevels of oppositely signed magnetic quantum number, the changes in the 
index of refraction for the two components have opposite sign. As a result of the index 
change, the components acquire a relative phase shift which leads to a rotation of the 
polarization direction which is proportional to the magnetic field for small fields. Thus, by 
measuring the rotation as a function of applied magnetic field, we obtain information about 
the dispersive properties of the medium. As the magnetic field increases, the absorption also 
increases because the splitting between the sublevels destroys the two-photon resonance, so 
the detection of the transmitted electromagnetic field intensity as a function of the magnetic 
field allows us to study the absorptive properties of the medium. Hence, using the nonlinear 
Faraday technique we can easily study both dispersive and absorptive properties of the EIT 
resonance simultaneously. 

We have measured the polarization rotation slope d<p/dB and the transmission I ou t/^in 
of the cell. The detailed shape depends on the particular cell and laser beam size, but a 
typical result is shown in Fig. ^. The individual points on this plot correspond to different 
atomic densities. These data cannot be fit by existing theoretical considerations for this 
system assuming a constant decay rate 70 ||15|| . This fact is demonstrated by extrapolating 



the low density end of the curve (where I out /Im is nearly 1) with a constant decay rate, as 
shown in the dashed curve and inset. However, taking into account the effect of radiation 
trapping we can understand this data quite well. 

We first analyze these data in terms of a simple theoretical model and then by a detailed 
numerical simulation. To simplify the analysis we neglect the process of optical pumping 
and assume a closed system. We include a dephasing of the ground-state coherence with rate 
7o and neglect population exchange between the ground states. As in previous treatments, 
we model the effect of radiation trapping by introducing an incoherent pumping rate R from 
ground states \b±) to excited state \a). Ris not a constant, but a function of all parameters 



of the system jy|. In particular, R depends on the excited-state population, because as 



more population is transferred to the excited state, more radiation will be produced that 
can eventually cause this incoherent excitation. We assume that R may be of the same order 
of magnitude as 70 but that it is much less then the radiative decay rate 7 r of transitions 
\a) -> |6±). 

We can understand the origin of this incoherent pumping by considering a two-level 
system coupled to a radiation reservoir. The reduced density matrix operator derived in 



Weisskopf-Wigner approximation ]16[ has the form 



p(t) = -n th lr [a-<r+p(t) - &+p(t)&J\ - 

(n t h + l)7r [(T+a-p(t) - a-p(t)a+] + H.c. , (1) 

where n t h is thermal average photon number in the reservoir, 7 r is the atomic decay rate of 
the upper level, er_ = \b)(a\ and cr+ = \a)(b\. To model atomic excitation by the incoherent 
radiation in the reservoir, the incoherent pumping rate can be written as R = 2^ r n t h • From 
Eq. (HD we find the equations of motion for the excited state population: 

p aa = -2j r (n t h + l)p aa + 2j r n th p bb . (2) 

In an optically thin atomic medium the probability of photon reabsorption is small and 
nth — 0- However, in optically thick media photons diffuse slowly and n t h 7^ 0. The value of 
nth can be estimated from a rate equation 

n th = -r e n th + r a p aa , (3) 

where r e is the photon escape rate and r a is the pumping rate due to the atomic decay. Both 
r e and r a depend on the geometry of the system and on the atomic density. In steady state, 
n t h = r a p aa /r e , and it follows from Eq. (g) and p aa + p bb = 1 that r e > r a . 

It is convenient to formally introduce a function f(N) defined by r a /r e = //(l + /) 
which characterizes the radiation trapping, such that f(N) > and f(N = 0) = 0. In 
the case when most population is in the ground state, p aa C 1, we see from Eq. (g) that 
Paa ~ ~ 27rPaa/(l + /)• I n the limit of low light intensity <C 7 r and weak radiation 



trapping n t h <C 1 the cw light propagation obeys d\Q\ 2 /dz ~ —2n^ T p aa , where Q is the Rabi 
frequency of the transition, k = (3/87r)iVA 2 7 r . , A is the wavelength, and z is the distance 
of propagation through the medium. We thus have a simple equation for the incoherent 
pumping rate due to radiation trapping: 

l^V)d 

k l + f(N) dz l 1 V ; 

This is a very intuitively appealing model: radiation trapping can exist only if the coherent 
radiation is absorbed by the system and is scattered due to spontaneous emission. 

We now return to the problem of radiation trapping in a three-level system with EIT. 
Because the Doppler distribution depends on atomic density (temperature), it is simplest to 
study radiation trapping in the Doppler- free limit of EIT, i.e., when the absorption and the 
dispersion do not depend on the width of the Doppler distribution Wd . Doppler averaging 



shows that this condition is fulfilled for relatively large light intensities \£l{z)\ 3> WdJ^o/'Yr 
for any z. 

The stationary propagation of the right and left circular polarized electric field com- 
ponents through the atomic vapor is described by Maxwell-Bloch equations in the slowly- 
varying amplitude and phase approximation ||16|| . We solve the equations by considering only 
the lowest order in 70, R and So, assuming |f2_(,2)| 2 « |fi + (,z)| 2 , where fl± are the complex 
Rabi-frequencies of the two optical fields. We separately consider the spatial evolution of the 
amplitudes and phases of these complex Rabi-frequencies by writing £l±(z) = |f2±(z)|e 1 ^ ± ( z - ) , 
and derive equations for the total intensity |f2| 2 and the relative phase = 0_ — <p + 

<l 'n\ 2 = -K(i + R), (5) 



dz 



To solve these equations we must specify the functional form of the incoherent pumping rate 
R in the three level configuration. From the general properties of radiation trapping || and 
from the results from radiation trapping in a two-level system, we assume that in the case 
of Doppler-free EIT the incoherent pumping can be modeled by Eq. (f|) as in the case of a 
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two-level system. As we shall see, this model works very well. 

With this form for R, Eq. (|5|) can be easily solved and we arrive at 



Q(z) 



2 



1 - loKZ (1 + f(N)) (7) 
|ft(0)| 2 



Q(0) 

so from Eq. (f|) we have R = f(N)j . Integration of Eq. for the phase yields 



d(p(z) 



dB 



2» B _ ln 



n(o) 



n(z) 



(8) 



Detection of |f2(L)/fi(0)| 2 and d<p(L)/dB, allows us to infer the value of the coherence decay 
rate as a function of the atomic density and estimate the radiation trapping effect. Thus we 
see that for optically thick media the coherence decay rate increases with the density. 



For smaller intensities \£l(z)\ <C Way^o/lr , Doppler-free EIT is not established, so the 
approximation Eq. @ is not valid and we do not discuss this regime here. 

Based on the low-density data in Fig. |2| (for which radiation trapping is negligible) 
Eqs.(0) and (|j) allow us to determine the coherence decay rate to be 7 ~ 0.0047 r . Given 
this value, we can then use the high density set of these data and Eq. (|j) to obtain the 
incoherent pumping rate R due to radiation trapping. The dependence is shown by the dots 
in Fig. |3|. In general, the functional form of R is not an "absolute" and it changes if the 
cell geometry or laser beam size changes, which is a key signature of the effect of radiation 
trapping. 

The probability of photon reabsorption becomes significant when the medium becomes 
optically thick on the length scale of the atomic cell size (under the Doppler-free EIT 
condition almost all atomic population is in the ground state), or 

± NX ^d^r>l. (9) 

8tt W d K ' 

For our experiment jr/Ap w 0.01, so Eq. (||) is fulfilled for IV > 5 x 10 10 cm -3 . For 
densities less than this, radiation trapping is negligible and we have R ~ 0. Above this 
value there are two distinct regimes of behavior, with both seen in Fig. [| At low density we 
have that R increases linearly with density due to photon absorption and emission within 



the cell. If the atomic beam is narrower than the radius of the atomic cell, as we have in 
our experiment (d ~ 0.1D), the next regime occurs for densities when photon reabsorption 
becomes significant inside the laser beam. In our case this is N > 5 x 10 11 cm -3 . 

To confirm our simple analytical calculations we have also made detailed numerical simu- 
lations of the experiment. We have considered light propagation in a thirteen- level Doppler- 
broadened system corresponding to the F = 2 — > F' = 1, 2 transition in 87 Rb. The decay 
of the atomic coherence was modeled by finite time of the flight through the laser beam 
(an open system). We solved the density matrix equations in steady state using the coher- 
ence decay rate as a fit parameter. In other words, we choose the effective coherence decay 
rate 70 + R in such a way that our numerical points for the dispersion dcj)/dB and intensity 

1 out I 'I'm corresponds to the experimental results. This is shown in the solid line in Fig. |2|. The 
dependence for R/^q obtained this way is shown in the solid line in Fig. |[ We see that the 
simple analytical analysis of the data coincides with the simulations for low atomic densities 
and diverges slightly for high densities. We explain this difference by inadequate intensity 
of the laser light. The maximum intensity of our laser (~ 100 mW/cm 2 ) corresponds to a 
Rabi frequency |Oq| ~ 3.67,., which lies on the edge of Doppler-free region determined by 
l^o I > Wdyjo/jr ~ 67,.. The absorption further decreases the intensity resulting in the 
Doppler broadening becoming important, unlike in our simplified calculations. 

Finally, we note that the observations reported here cannot be explained by spin exchange 
collisions between the atoms. The collisional cross section for Rb atoms is approximately 

2 x 10~ 14 cm 2 |10[ which results in a coherence decay rate 70 ~ 2 x 10 _5 7 r for the densities 
reported here. This is approximately two orders of magnitude less than the time-of-flight 
limited coherence decay rate 70 ~ 4 x 10 _3 7 r that we measured. 

In conclusion, we have shown both experimentally and theoretically that the effect of 
radiation trapping enhances the decay rate of the atomic coherence established by linearly 
polarized laser radiation between Zeeman sublevels. This effect leads to significant increase 
of the residual absorption in EIT experiments with optically thick atomic vapors. 
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FIGURES 

FIG. 1. Diagram showing the experimental setup; inset: Idealized three-level A-scheme con- 
sidered in the theoretical calculations. 



FIG. 2. The dependence of rotation rate d(j)/dB on transmission through the system I on t/hn '■ 
experimental (dots), previous theory with 70 = 0.0047 r (dashed line), and obtained by numerical 
simulation including radiation trapping, (solid line). 

FIG. 3. The incoherent pumping rate R/jo due to radiation trapping as a function of atomic 
density N: calculated by applying Eq. (||) to the data (dots) and obtained by numerical simulation 
(solid line). 
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Radiation trapping in coherent media 
Matsko et al. — Fig. 1 
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Radiation trapping in coherent media 
Matsko et al. — Fig. 2 
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Radiation trapping in coherent media 
Matsko et al. — Fig. 3 
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